Harmonic generation from YBa2Cu307_5 microwave resonators 
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Measurement of harmonic power on a suspended stripline microwave resonator patterned out of 
thin film YBa2Cu307_a (YBCO) has been carried out as a function of temperature and microwave 
power. The third harmonic power P3 shows a subtle nonlinear dependence on the fundamental power 
Pi on a log-log scale. Fit to a straight line yields slopes of P3 vs. Pi between 1.5 to 3, contrary 
to elementary calculations that predict a slope of exactly 3. It is shown that third harmonic power 
generated from higher order terms in the nonlinear impedance could account for the discrepancy 
quantitatively without resorting to any additional ad-hoc assumptions. Results of measurement of 
fifth harmonic power is presented as a confirmation. A small second harmonic power is also observed 
that does not show anomalous power dependence within the experimental accuracy. 



Since the day of its discovery, the new breed of ox- 
ide superconductors with transition temperatures above 
the boiling point of nitrogen have shown great potential 
as a replacement for metals in passive microwave struc- 
tures owing to their very low surface resistance 
However, widespread applications have been hampered 
by the inherent nonlinearity that these materials display. 
The nonlinearity manifests itself as an increase of inser- 
tion loss of the device with power level as well as har- 
monic generation and intermodulation distortion. A sys- 
tematic understanding of the phenomenology and micro- 
scopic origin of the nonlinear response is therefore very 
important for successful device applications. 

Measurement of the third order intercept (TOI), de- 
fined as the input power at which the fundamental and 
the third harmonic output power are equal to each other, 
is a quantitative measure of the nonlinear response of the 
device. The TOI can be measured both via a single tone 
(harmonic generation) as well as multi-tone (intermod- 
ulation distortion) measurement. Usually the former is 
carried out in a nonresonant transmission line structure 
and the later in a high-Q resonant structure with the ap- 
plied tones close to the resonance frequency /o- While the 
single tone experiment has the simplicity of design (only 
one signal source is needed), it lacks the sensitivity of 
the multi-tone measurement. The multi-tone measure- 
ment has the drawback of additional complexity owing 
to the need of multiple signal sources, power combiners 
and carefully matched waveguides. In the experiment 
described in this paper, we have carried out single tone 
measurements on a resonant structure, combining the 
advantage of both techniques. The problem of single- 
tone measurement on a resonant structure, namely, the 
shift of resonant frequency due to any change in exper- 



imental parameters was solved by dynamically tuning 
the synthesized source before taking each individual data 
point. Since A/ = h{P m i n ) - fo(P) « AX S where 
Xs = 27r/io/oA is the surface reactance, the effect of mi- 
crowave power induced change to the kinetic inductance 
of the device is also measured simultaneously. 

The experiments were carried out on films patterned 
out of five different YBCO wafers on LaA103 substrate 
obtained commercially. The wafers were patterned using 
standard photolithography techniques (wet etching) into 
straightlines 100/ito wide and 0.4/xm thick. Using a stan- 
dard microwave formula, the characteristic imedance Zq 
was evaluated to be 72 il. The patterned structure has 
YBCO only on one side, i.e. it does not use supercon- 
ducting ground planes. The copper package housing the 
film also provides the necessary grounding. This gives 
us the necessary sensitivity to measure the effect of very 
small applied DC magnetic fields which would be shielded 
out by a superconducting ground plane. Such a capabil- 
ity is important and unique to our setup as small DC 
fields have been shown to improve device characteristics 
[Q. The results of the field measurements, kinetic in- 
ductance effects, details of temperature dependence and 
some related topics will be described elsewhere ||. 

The resonators (fo — 4.3 GHz) were characterized us- 
ing a vector network analyzer. The unloaded Q was 
~ 1.2-1.8 x 10 4 at 10K. Coupling was subsequently in- 
creased to transmit maximum power through the res- 
onator. The nonlinearity manifests itself as a distortion 
of the Lorentzian shape of the resonance curve [^| and a 
decrease of Q. 

The synthesizer was set to generate a CW signal at the 
resonance frequency that was amplified by 30 dB, filtered 



*e-mail: srinivas@neu.edu, Web: http://sagar.physics.neu.edu/ 



1 



by a coaxial 5 GHz low pass filter and was fed to the res- 
onator. The output of the resonator was analyzed with 
a spectrum analyzer. 

The aspect of harmonic generation or intermodulation 
distortion that is most commonly investigated is its de- 
pendence on the input power. The knowledge of this de- 
pendence is very important both from the point of view 
of device design as well as understanding the mechanism 
of nonlinearity. A simple algebraic argument, as shown 
below, predicts the third order nonlinear power to scale 
as the cube of the input power. However, there is a lack 
of consensus about this in the experiments carried out 
by different investigators. Not only does the dependence 
differs from a cube @,[llj , it also depends upon the input 
power. It is often seen that the dependence is closer to 
P 3uJ oc if rather than if §. 

In this letter, we show that this discrepancy comes 
about because of: 

• Failure to take into account the power dependence 
of the insertion loss of the device. 

• The assumption that the entire third harmonic 
power comes from the cubic term d 2 Zs/di 2 of the 
current-dependent nonlinear impedance Z$- 

The first point, often neglected in the analysis of non- 
linear frequency conversion data is that the harmonic 
power is generated by the current flowing in the device 
which does not scale quadratically with the input power 
because of the power dependence of the insertion loss. In 
the case of a moderately nonlinear device this manifests 
as a distortion of the Lorentzian shape of the resonance 
in the frequency domain corresponding to the power be- 
ing transfered to higher harmonics. In the extreme case, 
this can lead to saturation and discontinuous changes in 
output power. Higher levels of applied power go to heat- 
ing the device decreasing the output power. Fig. ^ shows 
the value of the transmission amplitude S\2 for one of 
the resonators obtained from a vector network analyzer 
as the frequency is swept about the resonance at input 
power levels of 10, 5, 0, —5 and — 15dBm respectively, 
lOdBm corresponding to the lowest trace and — 15dBm 
corresponding to the highest. Decreasing the input power 
below — 15dBm value did not increase S\2 any further, 
indicating that the device is operating in the fully lin- 
ear regime. Similarly, increasing the input power above 
10 dBm did not seem to lower the saturation any futhcr. 
Depending on the sample, other manifestations of non- 
linear behavior, such as a step-like change superimposed 
on a Lorentzian S12 trace were also observed. Many of 
these can be modeled as a lumped resonator involving 
nonlinear elements |^|. Extremes of such nonlinear be- 
havior also serve as a characterization of the quality of 
the superconducting film as more granular films tend to 
demonstrate this behavior. 
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FIG. 1. Transmitted power vs. frequency about the reso- 
nance frequency at input powers of 10, 5, 0, -5 and -15 dBm, 
from bottom to top. 

Before each measurement of harmonic power reported 
below, the resonance was characterized using a network 
analyzer to ensure that the input/output antennas were 
not overcoupled and to verify that the response of the 
resonator was still the same. Direct coupling between 
the antennas were found to be at least 50 dB below the 
resonance, ensuring that the transmitted power comes 
from microwave currents induced in the resonator only. 

It is therefore clear that trying to calculate the slope 
of the harmonic power against the input power can give 
misleading results. One way to correct for the error is to 
calculate the RMS current in the resonator taking into 
account the insertion loss of the device (h]] . In this paper 
we plot the harmonic power against the output power at 
the fundamental frequency which automatically compen- 
sates for insertion losses. 

The second implicit assumption made in the analysis 
of the dependence of the third harmonic power on the 
input power is that the entire contribution to the former 
comes from the cubic term in the nonlinear impedance. 
However, any higher order odd term will also contribute 
to the third harmonic. Assuming v = ir\ + i 3 r 3 + i 5 rs, 
i = io sin Lot the third harmonic term can be written as 

V3w = 4* Zq (l + |^of|)i the negative sign implying 
that the 3rd harmonic is phase shifted by it from the 
fundamental. Since spectrum analyzer measurements are 
not sensitive to relative phase, we ignore the negative 
sign. Choosing a reference value for current i r (for ex- 
ample, the current that generates 1 mW of power in case 

of a dBm scale) and taking logarithm of both sides and 
■2 

assuming «C 1, we get 
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The presence of a fifth (or higher) order contribution to 
the third harmonic power would be indicated by a simul- 
taneous observation of the fifth harmonic signal. 
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Typical data for second, third and fifth harmonic of 
one of the films at 30 K is shown in fig.|[ The inset 
shows second and third harmonic data of a different film 
from a different manufacturer at 67.5 K. The second har- 
monic power is much weaker than the third (and even 
the fifth) as expected of a superconductor. The fifth 
harmonic power falls below the noise floor of the spec- 
trum analyzer at a fundamental power of ~ —8 dBm 
because the noise floor of the spectrum analyzer is sig- 
nificantly higher at higher frequencies. The second har- 
monic power, being very weak, is difficult to measure and 
analyze reliably. In the cases where a reasonably strong 
second harmonic signal was seen, its slope against the 
fundamental was close to 2:1 on a log-log scale, as is ex- 
pected. However, the third harmonic power has a much 
more complex dependence on the input power. Fitting 
the plot to a straightline yields a slope of anywhere be- 
tween 1.5 to 3 depending on the film and the temperature 
with a variance of a few %. Usually such values of the 
variance signifies that the fit is very good. However, as 
fig|| shows, the data has a richer structure than a straigh- 
line. 
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FIG. 2. Typical second, third and fifth harmonic power 
data at 30 K. Inset: Second and third harmonic power of a 
different film at 67.5 K 

It is therefore obvious that the third harmonic power 
cannot be coming from the third order nonlinear term 
only, as has traditionally been assumed. We can fit the 
data quantitatively by keeping one to two more term in 
the Taylor series expansion of the nonlinear impedance, 
depending on the sample. Fig.|| shows the data of fig.^ 
along with a straightline fit (with a slope of 2.3118 and 
variance of 0.0115) and a fit calculated by keeping the 
third harmonic contribution from a fifth and a seventh 
order term. The second fit is nearly indistinguishable 
from data. The top inset shows third harmonic data of a 
different film at 55 K, a straightline fit (with slope of 2.48 
and variance of 0.1142) and a fit keeping upto seventh or- 
der terms. The lower inset shows the second harmonic 



data of the same film at the same temperature as the up- 
per inset and a straightline fit with a slope of 2.001 and 
variance of 0.04). In those films (such as the one in fig.||) 
where a seventh order term was necessary to obtain a sat- 
isfactory fit, one should be able to see a seventh harmonic 
signal as well. We were, however, unable to measure it as 
it was outside the bandwidth of our spectrum analyzer. 
A similar method of analysis has also been used by other 
investigators H]. However, to our knowledge, this is the 
first measurement of the fifth harmonic which validates 
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FIG. 3. The same data as in fig.^. A straightline fit and 
one keeping a 5th and 7th order contribution is also shown. 
Inset: (Top) Third and (Bottom) Second harmonic data for 
a different film at 55 K along with fits. 

In principle the above type of analysis should be ap- 
plicable to the second harmonic power also. However, 
since the measured second harmonic data fits reasonably 
well to the predicted 2:1 behavior, it is likely that the 
contribution to the second harmonic power from higher 
harmonics is small. In any case, our low pass filter had 
a harmonic passband that coincided with the fourth har- 
monic frequency. We were, therefore, unable to measure 
this quantity. Let us point out here that the even har- 
monics do not contribute to the third harmonic power 
and odd harmonics do not contribute to the second har- 
monics. Physically speaking, this implies that time re- 
versal symmetry breaking and preserving mechanisms of 
harmonic generation are independent of each other. 

There are a few mechanisms for the observed second 
harmonic power, e.g. Josephson vortices, intrinsic mech- 
anisms breaking time reversal symmetry or magnetic ef- 
fects superimposed on the superconducting state. Har- 
monic generation experiments, being macroscopic in na- 
ture, cannot distinguish between these mechanisms and 
it is quite likely than several different mechanisms con- 
tribute to even harmonics simultaneously. 
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